Perpendicular magnetic anisotropy in nanoclusters caused by surface segregation and 

shape 



Stefan Heinrichs 1 , Wolfgang Dieterich 1 and Philipp Maass 2 
1 Fachbereich Physik, Universitat Konstanz, 784-57 Konstanz, Germany 
2 Institut fur Physik, Technische Universitat Ilmenau, 98684 Ilraenau, Germany 
(Dated: September 13, 2005) 

The growth of binary alloy clusters on a weakly interacting substrate through codeposition of two 
atomic species is studied by kinetic Monte Carlo simulation. Our model describes salient features 
of CoPt3-nanoclusters, as obtained recently by the molecular-beam epitaxy technique. The clusters 
display perpendicular magnetic anisotropy (PMA) in a temperature window of growth favorable for 
applications. This temperature window is found to arise from the interplay of Pt surface segregation 
and aspect ratio for cluster shapes. Conclusions are drawn how to optimize growth parameters with 
respect to PMA. 
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Growth of frozen-in metastable structures by molec- 
ular beam epitaxy (MBE) opens a wide field for pro- 
ducing materials with novel properties. In most cases, 
the growth conditions require a description based on far- 
from-equilibrium kinetics. Especially interesting are mul- 
ticomponent systems or systems with internal degrees of 
freedom (magnetic or electrical moments) that couple to 
the atomic structure. A rich variety of material proper- 
ties is then expected to emerge from the elementary pro- 
cesses of deposition, adatom diffusion, surface nucleation, 
as well as from structural and compositional changes of 
islands, from surface segregation or magnetic ordering. 

MBE growth of binary magnetic alloy thin films and 
nanoclusters has recently attracted much attention in the 
search for materials with an easy magnetic axis perpen- 
dicular to the substrate 0, 0] • Such perpendicular mag- 
netic anisotropy (PMA) can be utilized to increase mag- 
netic storage densities. In thin magnetic films the dipolar 
form anisotropy favors alignment of the magnetic mo- 
ments parallel to the substrate. To achieve PMA in thin 
films, the dipolar form anisotropy has to be overcome by 
effects based on local atomic arrangements. 

The system Co/Pt is particularly intriguing, as PMA 
could be obtained in different ways: (i) In Co/Pt multi- 
layers the magnetic moments of Co-atoms preferentially 
align along the interfacial Co-Pt bonds and PMA is found 
in magnetic X-ray circular dichroism studies to originate 
at the interfaces [3j . (H) In disordered films with CoPt3- 
stoichiometry grown in [111] direction, PMA is found in 
a temperature window 500 K < T < 820 K Q. XAFS 
and EXAFS experiments show the preferential oc- 
currence of Co/Pt bonds out of plane. These were as- 
sociated with surface-induced formation of flat Co-rich 
domains, as suggested by the increase of Curie temper- 
atures with respect to those for disordered bulk Co/Pt- 
alloys 0]. (Hi) In CoPt3 nanoclusters with 300 to 1000 
atoms, grown by codeposition of Co and Pt on WSe2 
van-der-Waals surfaces, PMA is observed near ambient 
growth temperatures pfl, favorable for applications. The 



underlying structural anisotropy, however, is not under- 
stood so far. 

In fact, the occurrence of PMA in nanoclusters is sur- 
prising in view of the kinetic limitations induced by com- 
paratively low growth temperatures. Several mechanisms 
can be envisaged in principle. One possibility might be 
a bulk anisotropy arising from formation of Co-rich do- 
mains in the bulk, similar to Q . Alternatively, PMA con- 
nected with the surface could be induced by Pt surface 
segregation. A further question concerns the interplay 
between LI2 ordering and PMA. 

In this Letter we investigate non-equilibrium struc- 
tures that emerge from the combined effects of depo- 
sition, surface diffusion and ordering phenomena in a 
model for binary alloy growth. We find that surface seg- 
regation together with cluster shape formation leads to a 
structural anisotropy that induces PMA. There is no evi- 
dence for significant bulk contributions facilitating PMA. 
In order to account for the experimentally observed phe- 
nomena, exchange processes during growth between low- 
coordinated atoms at the cluster surface have to be in- 
cluded. 

To demonstrate these mechanisms, we analyze a model 
based on nearest neighbor interactions Vaa, Vab, Vbb 
between atomic species A and B, which are adjusted 
to the CoPt3 system (A=Co, B=Pt). The parameter 
J = (Vaa + Vbb — 2Vab)/4 fixes the bulk order-disorder 
temperature T = l.83J/k B 8\. With T a = 958 K 
we take fcsTo/1.83 = 45 meV as our energy unit. The 
parameter h = Vbb — Vaa controls the surface segrega- 
tion of B atoms at equilibrium. Segregation with nearly 
100% Pt at the surface ^3 is compatible with h ~ 4. 
The parameter Vq — (Vab + Vbb)/2 represents the aver- 
age bond energy in the Ll2-ordered crystal, and together 
with J and h uniquely determines the bond energies Vaa, 
Vab > and Vbb ■ Its estimation from the cohesive energy of 
LI2 ordered CoPt3, however, would not be appropriate, 
because in real materials the total binding energy is not 
given by a sum of independent bond interactions |llj |. 
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FIG. 1: (color online) Structural anisotropy P and LI2 order 
parameter ^ for clusters of size 1000 for different substrate 
temperatures, with h — 4. The solid (dashed) lines are a 
guide to the eye for exchange barriers (7 X = (U x = 00). The 
inset shows a typical configuration ([/ x = 0, T = 1) before 
subsequent equilibration. Pt atoms are marked in blue (dark) 
and Co atoms in red. 



Rather we note that during cluster growth the mobility 
is effectively restricted to the surface, so that the number 
of relevant bonds per atom is between 3 (for a diffusing 
adatom) and 7 (for an atom at an edge). Parameters for 
a variety of corresponding processes were calculated for 
Pt with DFT 12] and yield V ~ -5, which we adopt 
here. 

Binding to the substrate surface is represented by a 
weak potential V s = — 5 restricted to adatoms in the first 
layer. Compared to a Pt(lll) surface with three bonds 
of typical strength Vo, this amounts to about 1/3 of the 
energy of a single Pt-Pt bond. 

To model the growth of the clusters, we perform kinetic 
Monte-Carlo (KMC) simulations on an fee lattice with 
periodic boundary conditions in the lateral directions. 
Compared to previous growth simulations of binary al- 
loys with solid-on-solid models or immediate downward 
funncling, our program is fully three-dimensional, so 
that for instance diffusion on equivalent side facets is 
automatically equivalent to diffusion on the top ter- 
race. A rejectionless continuous-time algorithm is imple- 
mented that generates cluster configurations with statis- 
tical weights according to the underlying master equa- 
tion. The possible elementary processes are jumps of A 
or B atoms to vacant nearest neighbor sites and the de- 
position of an atom. The jump of an atom occurs with a 
rate vexp(-U/kBT)mbx(l,exp[-AE/kBT]), where AE 
is the energy difference before and after the jump, U is 
an energy barrier and v an attempt frequency. Diffusion 
measurements of Pt on Pt(lll) give v ~ 5 x 10 12 s _1 
and U ~ 5. These values are taken to be the same for 
all jump processes. Deposition of atoms occurs with a 
rate F = 21 monolayers per second. The overall proce- 
dure automatically entails adatom diffusion, deposition 



and re-evaporation, nucleation, step edge diffusion, and 
step edge barriers. 

A typical cluster with N = 1000 atoms is shown in the 
inset of Fig.^ As in the experiments [7J, side facets are 
of {111} and {100} type and the aspect ratio between 
lateral and height dimensions is about 3:1. 

To quantify the long-range chemical ordering we de- 
termine the average intensity / of LI2 superstructure 
peaks and normalize it with respect to surrogate struc- 
tures, where the cluster shape is maintained, but the 
atomic configuration is changed. By occupying precisely 
one of the four sc-sublattices with Co atoms we ob- 
tain the intensity /id for perfect LI2 ordering, while a 
random occupation with the stoichiometric concentra- 
tion yields / ran . The order parameter is then defined 
as ^ = (I — Iran) /(/id — / ra n)- As can be seen from 
Fig. HI (triangles, dashed line), 4" decreases by lowering 
the temperature due to kinetic suppression of the equi- 
librium state, in accordance with experiments. The ex- 
trapolated onset temperature for Ll2-ordering, however, 
is much lower than in experiments. On the other hand, 
since the deposition rate in our simulations is larger than 
the experimental one, one expects a shift of calculated 
relative to experimental data to higher T. 

This problem can be resolved by allowing for direct ex- 
change processes between unlike neighboring atoms near 
the surface. As shown below, such processes also drive 
the Pt surface segregation which is crucial for producing 
PMA. In the model described before, only indirect ex- 
change processes between two atoms can take place via 
intermediate states, as a result of two or more direct ex- 
changes with vacancies. These are, however, extremely 
rare since realization of the intermediate states requires 
bond breaking steps with high activation energies. Direct 
exchange processes are found for Co deposited on Pt(lll) 
over a wide temperature range 250-520 K 0, Q, 0|. 
They are especially frequent for low-coordinated atoms 
on top of terraces (with coordination 3) or at step edges 
(with coordinations 4 or 5). 

To include these processes in the model, we allow a pair 
of unlike neighboring atoms to exchange sites if one of 
them has a coordination in the range 3 to 5 and the other 
one in the range 8 to 10. In addition to the difference 
between initial and final energies, an increased exchange 
barrier U x + U has to be surmounted. Setting U x = as 
the opposite extreme to U K = 00 (no exchange processes), 
we see in Fig. 1 that the onset of Ll2-ordering occurs 
near T ~ 1 (corresponding to T — 523 K), which is in 
reasonable accordance with experiments. 

We now investigate the possible occurrence of a struc- 
tural anisotropy that can induce PMA. Regarding mag- 
netic properties in Co-Pt alloys, we adopt the view- 
point, suggested by the experiments discussed above 
3, 4], that the crystalline magnetic anisotropy origi- 
nates mainly from Co-Pt bonds, which tend to align the 
Co moments along the bond. PMA will therefore be 
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traced back to a prevalence of out-of-plane relative to 
in-plane Co-Pt bonds. Denoting the corresponding to- 
tal numbers of bonds in a cluster of N atoms by n2° pt 
and n^ oPt , the pertinent anisotropy is introduced as 
P = 05 oPt - n^ oPt )/N. PMA requires that P > 0. 

In the absence of exchange processes (U x = oo), sim- 
ulated values for P turn out, however, to be negative as 
shown in Fig. 1 . The reason is that the clusters are rela- 
tively flat and that surface segregation of Pt is weak. As 
a consequence, in a clusters' surface there exist more Co- 
Pt bonds in-plane than out-of-plane, while atoms in the 
inner part of the cluster give no contribution to P. The 
flat cluster shape, on the other hand, would favor PMA 
if strong Pt segregation, as realized at equilibrium, could 
build up. Indeed, when including exchange processes, the 
Pt segregation gets strongly enhanced and P can become 
positive. Actually, the cluster displayed in the inset of 
Fig. 1, showing pronounced Pt surface segregation, was 
computed with U x = 0. As seen from the figure, for 
U x = 0, P changes sign near T = 0.58 and reaches a 
maximum at an optimum temperature T max ~ 0.8. As 
in experiments 0] , the subsequent decrease of P is accom- 
panied by the onset and subsequent rise of Ll2-ordering. 

To investigate the mechanisms of PMA in more detail, 
we show in Fig.Et P as a function of cluster size N during 
growth for 4 different temperatures, and for a more real- 
istic nonzero exchange barrier U x > 0, which we set equal 
to the adatom diffusion barrier, U x = U = 5. Except for 
the two lowest temperatures, where Pt segregation is not 
sufficiently strong, P indeed becomes positive except for 
very small N (not shown) . Again there exists an optimal 
temperature T max ~ 1, where PMA is expected to be 
strongest. Changes of P during subsequent equilibration 
under zero flux are less then 15 % of the values near T max 
and will be ignored in the following. 

The existence of an optimal temperature, however, is 
not directly caused by LI2 ordering as one might conjec- 
ture from Fig. 1. Rather, it results from two competing 
effects: As shown in Fig. increasing T facilitates Pt 
segregation, leading to larger P. On the other hand, in- 
creasing T also drives the cluster morphology closer to an 
equilibrium with less flat shapes with smaller P. This ef- 
fect is shown in Fig. [2k, where we plot the T-dependence 
of the aspect ratio, defined by the ratio l\\/l± of the gy- 
ration radii of cluster sizes in the direction parallel and 
perpendicular to the substrate. While the aspect ratio 
does not depend on U x , the Pt segregation is significantly 
enhanced for lower U x , see Fig.^D. These results can ex- 
plain why there exists a temperature window for PMA 
in experiments and why in films with tall rather than 
flat clusters, parallel instead of perpendicular magnetic 
anisotropy is favored (lti} . 

Since the positive values of P arise due to the surface 
segregation of Pt, we expect n ( l oPt — n^ oPt ~ A 2 / 3 , and 

accordingly P ~ A -1 / 3 . As shown by the solid lines in 
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FIG. 2: (a) Anisotropy parameter P averaged over 20 clusters 
versus size N at different temperatures, with (7 X = 5. For the 
two highest temperatures, the proportionality P oc iV -1 / 3 is 
well established, see the lines, (b) Pt concentration Cpt in the 
outer shell of the cluster for various U x (symbol assignment as 
in (cj). (c) Ratio lu/lx of cluster gyration radii as a function 
of temperature. The line is drawn as a guide to the eye. 

Fig. this behavior is well obeyed by the data sets for 
T = 1 and T = 1.2. 

We are now in the position to estimate the magni- 
tude of the magnetic anisotropy E s = E{fl in plane} - 
E{jl out of plane}, as expected from the data in Fig. |2] 
E{p\ denotes the anisotropic part of magnetic interac- 
tions as a function of Co moments fl for saturated magne- 
tization. Within a simple bond picture we suppose that 
each Co-atom experiences a local crystalline anisotropy, 
which is a sum over the corresponding twelve Co-a bonds 
in the fcc-lattice with a — Co, Pt or V(vacancy), speci- 
fying the nearest neighbor occupation. Each bond with 
bond-vector S contributes a term — A a (fl ■ S/\fl\\5\) 2 to 
E{jl}, with anisotropy constants A a to be deduced from 
experiment. Note that for saturated magnetization a 
nearest-neighbor contribution to the anisotropic part of 
dipole-dipole interactions has exactly the same form, so 
that these are included already in the coefficients A a . As 
discussed before, the most important contribution arises 
from Co-Pt bonds. Anisotropy constants for Co-Pt in 
multi-layer systems and wedges were measured for in- 
stance with magnetic torque or the magneto-optical Kerr 
effect From these experiments we deduce A CoPt = 
160 fieV as a representative value. Similarly, from mea- 
surements of Co- vacuum interfaces and theoretical inves- 
tigations of a freely standing Co-monolayer 01 we es- 
timate A CoV ~ —67 fieV. The remaining parameter is 
A CoCo for which we retain only the dipolar contribution. 
Using n Co = 1.7/Xb El and \8\ = 2.72 A in CoPt 3 , 
this yields ,4 CoCo = 23 fieV. 

Within this bond model it is straightforward to cal- 
culate the anisotropy energy E s , which can be written 
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FIG. 3: (a) Total magnetic anisotropy energy E tot for 
N = 10 3 in dependence of temperature for different exchange 
barriers £/ x for exchange processes, (b) Surface and bulk con- 
tributions E B , -Edip in dependence of cluster size N. The cor- 
responding straight lines in the double logarithmic plot have 
slopes 2/3 and 1. Data for E to t are fitted according to Eq. @. 



how PMA can result from the combined effect of surface 
segregation and flat cluster morphologies. The interplay 
between segregation and aspect ratio of cluster shapes ex- 
plains the occurrence of a temperature window for PMA. 
The effect discussed here primarily relies on surface con- 
tributions and favors PMA, in contrast to the dipolar 
form anisotropy. Stronger segregation (e.g. due to lower 
exchange barriers) as well as smaller adsorbate-substrate 
and adsorbate-vacuum surface tensions (leading to flat- 
ter shapes) should broaden the temperature window and 
may guide the search for optimized materials. We expect 
that for T ~ T max there exists an optimal cluster size 
N opt for PMA and a critical cluster size iV c = (3/2) 3 iV p t , 
above which PMA disappears. It should be interesting 
to test these predictions in experiments. 

We thank M. Albrecht, M. Maret, J. Schafer and 
G. Schatz for helpful discussions and the Deutsche 
Forschungsgemeinschaft (SFB 513) for financial support. 



as , 

Es = t £ (A CoQ -^ CoV )(n5 OQ -nf oV ) 

a=Pt,Co (1) 

= K S N 2/3 

The last equality in QJ holds on account of P cx iV -1 / 3 , 
see Figs. 2a and 3b. Similarly, we found that (n5° Co — 
n^° c °) cx 7V 2/3 , but this term related to Co-Co bonds 
turns out to be significantly smaller than the first term 
in in view of our estimates for A Coa . Altogether we 
find K s ~ 200 /ieV at T = 1 and U K = 5. 

So far we neglected contributions to the anisotropy en- 
ergy from dipolar interactions between non-neighboring 
moments. These interactions favor an orientation of the 
magnetization parallel to the (lll)-plane and hence give 
a negative bulk contribution E d i p = —Kdi p N to the total 
anisotropy energy 

E tot = K S N 2 / 3 - K dip N (2) 

From calculated cluster structures, Co magnetic mo- 
ments 1.7/ib and induced Pt magnetic moments 0.3/ie 
HEH we obtain K dip -5.6 fieV at T = 1. This yields 
an jEtot that is displayed in Fig. 3a as a function of T 
for different J7 X and N = 1000. An optimal temperature, 
where .Etot is maximum, can clearly be identified. For 
example, for U x = 5, T max ~ 1. 

In Fig. we show E to t as a function of N for U x = 5 
and T = 1 ~ T max , together with the surface and dipolar- 
contributions E s cx TV 2 / 3 and E d [ p cx N. The solid line 
for Etot represents Eq. J5J, which predicts an optimal 
cluster size iV pt ~ 13500 with maximum Etot- The PMA 
is expected to disappear for N > N c ~ 5 • 10 4 . 

In summary, by simulating the fully three-dimensional 
kinetics of binary alloy cluster growth, we have shown 
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